Anomalous suppression of the superfiuid density in the Cu a Bi 2 Se3 superconductor 

upon progressive Cu intercalation 
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CuzBi2Se3 was recently found to be likely the first example of a time-reversal-invariant topo- 
logical superconductor accompanied by helical Majorana fermions on the surface. Here we present 
that progressive Cu intercalation into this system introduces significant disorder and leads to an 
anomalous suppression of the superfiuid density which was obtained from the measurements of the 
lower critical field. At the same time, the transition temperature T c is only moderately suppressed, 
which agrees with a recent prediction for the impurity effect in this class of topological superconduc- 
tors bearing strong spin-orbit coupling. Those unusual disorder effects give support to the possible 
odd-parity pairing state in CuiBi2Se3. 

PACS numbers: 74.20.Mn; 74.20.Rp; 74.25. Bt; 74.62.Dh 



A topological superconductor (TSC) is a supercon- 
ducting analog of topological insulators (TIs) and is char- 
acterized by a nontrivial topological structure of the 
Hilbert space, which is specified by nontrivial Z or Z-2 
indices.— ~— Its hallmark signature is the appearance of 
surface Majorana fermions, which are their own antipar- 
ticles and are of fundamental intellectual interest 4 Re- 
cently, it was theoretically predicted 9 - and experimentally 
confirmed^ that a superconducting doped TI material 
Cua;Bi2Se3 is likely the first concrete example of a time- 
reversal-invariant TSC. Its parent TI material Bi2Sc3 
consists of basic crystallographic units of Se-Bi-Se-Bi-Se 
quintuple layers, which are weakly bonded by the van der 
Waals force. Upon intercalation of Cu into the van der 
Waals gap, superconductivity appears below the critical 
temperature T c of up to ~ 3.8 K^i This is high for its low 
charge carrier concentration of only ~ 10 20 cm~ 3 Ai As 
a "superconducting topological insulator" , this material 
has attracted a great deal of interest.— 

Unfortunately, Cua;Bi2Se3 has a materials problem in 
that samples with a large superconducting volume frac- 
tion are difficult to obtain with the usual melt-growth 
method ) 11 ' 12 which has hindered detailed studies of the 
superconducting properties of this material. However, 
this problem has been ameliorated recently by the devel- 
opment of an electrochemical synthesis technique^ which 
allowed the synthesis of superconducting samples with 
the shielding fraction exceeding 50% Such an improve- 
ment made it possible to perform point-contact spec- 
troscopy on a cleaved surface of Cu2:Bi2Se3 to find a sig- 
nature of the Andreev bound state in the form of a pro- 
nounced zero-bias conductance peakji - which gives evi- 
dence for unconventional superconductivit y 27 ' 28 Know- 
ing that the symmetry of this materia l 29 ' 30 allows only 
four types of superconducting gap functions^ and that 
all possible unconventional states are topological^ it 
was possible to conclude that Cua;Bi2Se3 is most likely a 
TSC±P. 

Although the point-contact spectroscopy elucidated 
the possible TSC nature of Cu.rBi2Se3, the electron mean 
free path £ in the superconducting samples of this ma- 



terial is comparable to the coherence length £or^ ac- 
cording to the common belief ; 31 ' 32 the odd-parity pairing 
should be strongly suppressed by impurity scattering in 
such a situationj 33 ' 34 In this context, a recent theory by 
Michaeli and Fu addressed this issued 3 , and showed that 
odd-parity superconductivity in strongly spin-orbit cou- 
pled semiconductors such as Cua;Bi2Se3 are much more 
robust against the pair-breaking effect induced by impu- 
rity scattering than in more ordinary odd-parity super- 
conductors. Therefore, thanks to the role of spin-orbit 
coupling, T c of Cua;Bi2Se3 is expected to be rather in- 
sensitive to nonmagnetic impurities, which is similar to 
conventional superconductors^ 

In this Rapid Communication, we address the issue 
of disorder effects in Cu x Bi2Se3. Through our system- 
atic studies of the effects of Cu intercalation in this sys- 
tem, it turned out that increasing the Cu content beyond 
x ~ 0.3 in the superconducting regime does not increase 
T c or the carrier concentration, but its main effect is to 
enhance the residual resistivity po- This suggests that 
one can consider the Cu content x to be a parameter 
to control the disorder while keeping other fundamental 
parameters virtually unchanged. By looking at the data 
from this perspective, the x dependence of the superfiuid 
density obtained from the lower critical field shows an 
unusual disorder dependence that is distinct from that 
in conventional BCS superconductors, which gives sup- 
port to unconventional pairing. In addition, we show 
that the x dependence of T c is essentially a reflection of 
the disorder effect and is consistent with the particular 
odd-parity paring state that is supposed to be realized in 
Cu K Bi 2 Se 3 . 

00,312863 single crystals of slab- like geometry with 
various Cu contents 0.11 < x < 0.50 were prepared by the 
electrochemical technique described carlierJ^ The typi- 
cal sample size was 4 x 2.5 x 0.3 mm 3 . The magnetic 
field dependence of the magnetization, M(B), were mea- 
sured with a commercial superconducting quantum inter- 
ference device (SQUID) magnetometer (Quantum Design 
MPMS) with particular attention being paid to the low 
field regime^ Roughly half of the samples were also char- 
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FIG. 1: (Color online) Cu content x dependences of (a) crit- 
ical temperature T c , (b) residual resistivity po, (c) supercon- 
ducting shielding fraction, and (d) normal-state carrier den- 
sity n. For comparison, the carrier density in pristine Bi2Se3 
(< 10 19 cm -3 ) is shown with a diamond. Note that n is essen- 



tially x independent and remains low at 
The dotted lines are guides to the eyes. 
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acterized by transport measurements by a standard six- 
probe method. Figure [T] summarizes the x dependences 
of T c , po (defined as p at T = 5 K), the superconducting 
shielding fraction at T = 1.8 K, and the charge carrier 
concentration n (determined from the Hall coefficient at 
5 K). Most notably, po strongly increases for x > 0.3 and 
n is basically independent of x at n ~ 1.5 x 10 20 cm" 3 i2 
Before presenting magnetic properties, we define and 
summarize important parameters. The layered structure 
of CuzE^Sea leads to anisotropics in the superconduct- 
ing parameters, and we denote the lower and upper crit- 
ical fields for magnetic fields parallel and perpendicular 
to the crystallographic ab planes as B c i )0 &, -B c i,c, B c2t ab, 
and -B C 2,e, respectively. Also, the penetration depths and 
the coherence lengths along the in-plane and out-of-planc 
directions are denoted as X ab , A c , £ a b, and £ c , respec- 
tively. The anisotropic Ginzburg-Landau parameters are 
defined as n ab = yj 'X ab X c / '(£ Q b£ c ) and k c = \ ab /£ ab M~— 
The upper critical fields are related to the coherence 
lengths via B c2 ^ ab = $o/27r£ a &£c and B c2 . c = $ /27r^ b 
with the flux quantum $o- In the Ginzburg-Landau 
theory, B c \ is related to the vortex line energy E via 
B c \ = Att^lqE/^o^ for extremely type-II superconduc- 
tors with k 1, one obtains E w [$1/(4tt\) 2 ] ln/c. How- 
ever, to take into account the vortex core energy, the 
In k term has to be corrected by adding 0.5^—, and the 
formula for B c \ ab becomes 



$0 1 

Bd,ab = —[ln(K a b) + 0.5]- — 

47T A a bA c 



(1) 



Hence, to calculate the Ginzburg-Landau parameter 
K ab , we use B cl>ab /B C 2,ab = {^K ab + 0.5)/2K 2 b . The 
anisotropy factor is defined as 7 = B c2 ^ ab / B c2 ^ c = 
B c i,c/B c i ab = X c /X ab and the penetration depths are de- 
termined by solving Eq.[T]for X ab by using A c = 7A Q (,. For 
the following discussion, we define the averaged penetra- 
tion depth A av = y A 2 b A c , which allows the calculation 



FIG. 2: (Color online) (a) Initial low-field M cxp vs B curves 
for a sample with x — 0.26 at various temperatures, (b) Re- 
duced magnetization AM after subtracting the initial linear 
Meissner contribution My m . The deviation points marked by 
arrows indicate B\ at each temperature, (c) Plot of B\ vs T 
for x = 0.26 together with a fit to the data, (d) Plots of AM 
vs 6 at 1.8 K for various x. 



of the superfluid density via n s = m*/(/ioe 2 A 2 v ) with the 
effective mass m* assumed to be x independent^ 

Figure 2 describes how B c i,ab is determined from the 
magnetization data M" eX p, which is essentially the same 
as was done in Ref. [l7l Figure 2(a) shows M exp (B) 
curves for a sample with x = 0.26 at various tempera- 
tures between 1.8 and 3.0 K; the T c of this sample was 
^3.2 K. The dashed line My m is a fit to the low- field 
magnetization at 1.8 K, representing the initial Meissner 
screening. Determining such a linear part for each tem- 
perature and subtracting it from the M exp (B) data yields 
AM = M ex p— Mii n which is shown in Fig. 2(b) (the data 
are shifted for clarity). The arrows mark the last field 
above which AM shows an obvious deviation from zero 
(i.e. AM becomes > 0.02), signaling the entry of vor- 
tices and defining the flux entry field B\ (T) . These data 
points arc plotted as B\ vs T in Fig. 2(c) and are fitted 
with the empirical formula Bi (T) = Bi(0)[l-(T/T c ) 4 ]4& 
Sometimes the B\{T) data scatter around the fitting line, 
which leads to a sample-dependent error bar on _Bi(0). 
To show the trend of how B\ changes with the Cu con- 
centration x, the 1.8-K magnetization data for various x 
are shown in Fig. 2(d); this plot clearly suggests that 
Bi(0) becomes systematically smaller for larger x. 

To determine B c \ from B±, one should consider the in- 
fluences of the demagnetization effect, the surface quality 
(Bean-Livingston surface barrier), and the bulk pinning 
effects. The surface barrier is not effective in rough sur- 
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FIG. 3: (Color online) x dependences of (a) B c i(0) for B || ab 
after the demagnetization correction, (b) corresponding pen- 
etration depth A av , (c) £> C 2(0) for fields parallel and perpen- 
dicular to the ab planes, and (d) anisotropy factor 7. Dotted 
lines are guides to the eyes. 



faces, which is the case in all of our as-grown samples 
(see Fig. S4 of Refj 1 ^) irrespective of the x values, and 
the bulk pinning is also extremely weak in CuzE^Sea as 
indicated by magnetic hysteresis data i 17 ' 47 As for the de- 
magnetization effect, albeit small (< 10%) in the present 
case, we have corrected for it by using the Brandt's for- 
mula for slab-shaped samples with an aspect ratio fr/a4£ 
B cl (0) = Si(0)/tanh ^/0.36 6/a. The obtained B cl for 
all samples are plotted vs x in Fig. 3(a). The correspond- 
ing A av values arc shown in Fig. 3(b); for calculating A av , 
we need the anisotropy factor 7 which was obtained from 
anisotropic B C 2 determined from the resistive transitions 
in magnetic fields applied parallel and perpendicular to 
the ab plane [Fig. 3(c)]. The obtained 7 is essentially in- 
dependent of x [Fig. 3(d)], which supports the idea that 
the main effect of Cu intercalation beyond x ~ 0.3 is to 
enhance the disorder without changing band structure or 
mobile carrier density. 

As already mentioned, the averaged penetration 
depth A av directly gives the supcrfluid density n s = 
m*/(tioe 2 A av ).— We normalize this value with the 
normal-state carrier density n, and Fig. 4(a) summa- 
rizes the x dependence of nf p = n s /n. One can see that 
n s /n is already only 0.3 at x ~ 0.10 where the supercon- 
ductivity starts to appear, and it is further suppressed 
with increasing x. This behavior is obviously a reflection 
of strong disorder caused by Cu intercalation that can 
be inferred in Fig. 1(b). Since it is known that disor- 
der causes a reduction in n s even in conventional BCS 
superconductors , 50 ! 51 it is prudent to discuss this behav- 
ior quantitatively. 

According to Anderson's theorem^ 5 - the superconduct- 
ing gap Ao and T c of conventional superconductors are 
relatively insensitive to small concentrations of nonmag- 
netic impurities. However, the supcrfluid density, which 
reflects the "rigidity" of the electronic system to electro- 
magnetic perturbations, is affected by disorder in con- 
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FIG. 4: (Color online) (a) Normalized superfluid density n s /n 
vs x, with n the normal-state carrier density. For symbols 
without error bars, the estimated errors are smaller than the 
symbol size, (b), (c) Semi-log plots of n s /n (b) and T c (c) vs 
the mean free path £. The upper axis gives the corresponding 
k^l value. To facilitate comparisons between (a) and (b) [(c)], 
data for different x values are indicated by different symbols. 
The solid line n? cs in (b) gives the expected disorder-induced 
suppression of n a /n for a conventional BCS superconductor. 
The inset shows nf cs -> 1 in the clean limit £ — > 00. The 
dotted lines in all panels are guides to the eyes. 



ventional superconductors.— ~— Indeed, the disorder de- 
pendence of n s has been studied in Nb and Pb and was 
found to follow the theoretical prediction ! 53 ' 54 We there- 
fore compare the disorder dependence of n s observed 
in Cu;,.Bi2Se3 to the expectation for ordinary BCS su- 
perconductors. For such a comparison, one needs to 
parametrize disorder, which is usually done by evaluat- 
ing where kp = y/3ir 2 n is the Fermi wave number— 
and t = hk-p / '(pone 2 ) is the mean free path.— 

For a pure BCS superconductor, the penetration depth 
in the 0-K limit is given by A^(0) = m*/(/ioe 2 n), be- 
cause n s is equal to n in the clean limit. In the presence 
of disorder, this Al(0) in the clean limit is modified to 
an effective penetration depth which is evaluated at T 
= K as Abcs(0) = A L (0)^1 + ^/^ > A L (0) in the 
local limit ; 57 ' 58 where £0 = Hvf/^Ao) is the Pippard co- 
herence length for pure superconductors (v-p = hkp/m* 
is the Fermi velocity and A is the BCS gap)^ From 
this Abcs(0) we calculate the superfluid density nf cs (£), 
which gives the disorder-induced suppression of n s for a 
conventional BCS superconductor. 

Figure 4(b) shows the comparison of the I dependences 
of rig Xp and nf cs (kp£ value is shown in the upper axis). 
In this figure, the BCS calculation is shown as a solid 
line and the inset shows the saturation of ??^ cs — > 1 in 
the clean limit I — > 00. Clearly, n| xp does not agree with 
, although both are suppressed with decreasing £, 
the suppression is much stronger in Cua;Bi2Se3 than is 
expected for a BCS superconductor. Also, it is useful 
to compare the result shown in Fig. 4(b) to that in Fig. 
1(b): At x > 0.3, the residual resistivity starts to increase 
drastically and I becomes shorter than ^25 nm; however, 
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ng Xp tends to saturate in this dirtier range of I < 25 nm. 
Moreover, for £ < 4 nm, n c s xp intersects the n^ cs curve. 
Hence, both the strong suppression in the intermediate 
disorder regime and the saturation tendency in the dirt- 
ier regime are anomalous. Such an anomalous behavior 
of n s /n is the main result of this work, and it naturally 
points to an unconventional nature of the superconduct- 
ing state in Cu^E^Sea. 

In contrast to the highly anomalous behavior of 
n° xp , the modest suppression of T c shown in Fig. 
4(c) resembles the behavior of dirty conventional 
superconductors^ One might hasten to conclude that 
such an ordinary disorder dependence of T c speaks 
against the odd-parity pairing, because the common be- 
lief for odd-parity superconductors is that T c is quickly 
suppressed with impurity-induced disorder ! 31 ' 32 How- 
ever, as we already mentioned above, the particular type 
of odd-parity pairing that is considered to be realized 
in Cu^E^Ses (Refs. I9lll0l) belies this common belief. 
This point was recently shown by Michaeli and Fu^ 3 - 
who analyzed the novel inter-orbital, odd parity stated 
proposed for Cu^E^Ses. The odd-parity pairing takes 
place between two p z orbitals with different parity at the 
upper and lower ends of the quintuple layers via attrac- 
tive short-range interactions. In such a state, the cru- 
cial disorder-induced pair breaking effect is significantly 
suppressed as a result of strong spin-momentum lock- 
ing. The dephasing rate of the Cooper pairs depends 
on the ratio of band mass and chemical potential, m//i; 
as this ratio becomes smaller, the superconductivity be- 
comes more robust. For Cuj;Bi2Sc3, this ratio has been 



estimated^ to be ~ 1/3 and the calculated T c depends 
only weakly on the impurity-induced disorder^ 3 - in quali- 
tative agreement with Fig. 4(c). Therefore, the observed 
disorder effect in T c is not inconsistent with the odd- 
parity pairing. 

To summarize, we report an anomalous suppression of 
the superfluid density n s /n probed by the lower critical 
field as a function of the Cu content x. Since it appears 
that the main effect of Cu intercalation beyond x ~ 0.3 is 
to enhance disorder without significantly changing band 
structure or carrier density, our result reveals the impact 
of disorder on the superconducting state in Cu^T^Ses. 
Most strikingly, in the intermediate range of disorder, 
n s /n is much more strongly suppressed than is expected 
for a dirty conventional BCS superconductor, while in 
the strongly disordered regime n s /n tends to saturate. 
In contrast, the occurrence of superconductivity itself is 
robust against disorder as indicated by an only moderate 
suppression of T c with x. The obviously anomalous be- 
havior in n s /n points to an unconventional pairing state, 
and the ostensibly normal behavior in T c is consistent 
with the theoretically-proposed odd-parity pairing state 
with strong spin-momentum locking. Altogether, our re- 
sult gives support to the possible odd-parity pairing state 
in Cu x Bi2Se3. 
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